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ABSTRACT 


Investigations  of  techniques  for  describing  and 
enhancing  the  accuracy  of  optical  linear  algebra  processors 
have  been  conducted.  Significant  accomplishments  include; 
(1)  development  and  simulation  of  a  system  model 
incorporating  device  dynamic  range  for  better  quantitative 
assessment  of  error-correction  code  performance;  (2) 
extension  of  earlier  statistical  models  to  include 
crosstalk,  background,  avalanche  gain,  flicker  and 
generation-recou±)ination  noise  effects;  (3)  construction  of 
the  Optical  Analysis  Simulation  Interactive  System  (OASIS) 
software  for  the  acquisition,  analysis  and  manipulation  of 
experimental  data  and  (4)  identification  of  major  noise 
sources  of  experimental  concern  using  OASIS. 
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RESEARCH  OBJECTIVES 


During  the  second  grant  funding  period  (March  1, 
1992-February  28,  1993)  ,  the  major  research  objectives 
have  been  to  continue  both  analytical  and  experimental 
investigations  on  techniques  for  enhancing  the  accuracy 
of  optical  linear  algebra  processors.  This  work  has 
been  divided  into  three  separate  but  related  areas:  (1) 
an  investigation  of  error  correction/detection 
techiniques  for  use  in  quantized  analog  optical  linear 
algebra  processor  applications  (2)  complete  statistical 
modeling  of  a  generic  optical  vector-matrix  multiplier 
taking  into  account  noise  effects  in  sources,  spatial 
light  modulators  (SLMs)  and  detectors  to  predict  overall 
noise  performance  and  (3)  experimental  characterization 
of  sources,  SLMs,  and  detectors  via  a  iow-noise 
measurement  facility  and  associated  data  acquisition  and 
analysis  software.  Details  of  these  investigations  are 
presented  in  the  following  sections  and  in  the 
publications  referenced. 


SUMMARY  OF  RESULTS 


Since  most  of  the  results  obtained  under  the  grant  are 
promptly  submitted  for  publication,  and  are  also  presented 
at  national  and  international  scientific  meetings,  we  will 
summarize  the  major  results  obtained  in  this  section,  with 
references  to  the  appropriate  journal  articles  and 
conference  proceedings, 

1 .  Encoding/Decoding  for  Error  Correction  in  Optical 
Computing 

This  phase  of  the  program  involves  finding  and/or 
developing  error-correction  codes  that  can  be  applied  to 
optical  linear  algebra  processors  (OLAPs)  to  enhance  their 
accuracy.  Error-correction  codes  have  been  successfully 
applied  to  increase  the  reliability  of  data  transmission 
systems.  The  codes  introduce  a  controlled  amount  of 
redundant  information  into  the  transmitted  data.  At  the 
receiver,  the  redundant  information  allows  errors  in  the 
received  data  to  be  detected  and/or  corrected,  thus 
improving  the  reliability  of  the  transmission  system.  In  a 
similar  manner,  this  idea  can  be  applied  to  OLAPs [1, 2, 3, i ] . 
Consider  an  optical  matrix-vector  multiplier  (OMVM) . 
Redundant  information  can  be  introduced  into  the  matrix  by 
increasing  the  rov.'  dimension  of  tne  matrix.  The  additional 
rows  in  the  matrix  are  formed  as  linear  combinations  of  the 
rows  of  the  original  matrix.  These  new  rows  provide  the 
redundant  information.  The  augmented  matrix  is  now  post 


multiplied  by  the  input  vector  to  form  the  product  vector. 
The  product  vector  contains  redundant  information  that  can 
be  used  to  correct  errors  in  the  elements  of  the  product 

vector.  Errors  can  occur  during  the  matrix- vector 
multiplication  due  to  noise  in  the  system  or  due  to 

component  failures. 

A  block  diagram  of  an  OMVM  employing  an  error- 
correction  code  is  shown  in  Figure  1.  In  general,  an  OMV.M  is 

designed  to  implement  the  matrix- vector  multiplication 

Y  =  Ax 

(1) 

where  A  is  an  M  X  N  matrix  and  x  is  an  N  X  1  vector. 

Redundancy  is  introduced  into  the  matrix  by  premultiplying 

the  matrix  by  the  generator  matrix  G  for  an  error  code,  or 
Aq  =  GA.  The  encoder  performs  this  operation.  The  generator 
matrix  G  is  an  X  M  matrix 


OMVM 


Figurf^  1  .  Block  diagram  of  an  optical  matrix-vector 
multiplier  employing  an  error-correction 
code . 


with  >  M.  Next  the  product  A^x  =  GAx  =  r’  is  computed  in 
the  OMVM.  The  product  r'  may  contain  errors  as  discussed 
earlier.  Note  that  the  multiplication  GAx  =  Gy  is  equivalent 
to  encoding  the  vector  y  as  is  done  for  data  transmission. 
Therefore,  the  same  decoding  algorithm  is  applied  to  r'  as 
for  the  case  of  data  transmission .  The  decoder  performs  the 
error  detection  and  correction  operations  to  produce  y.  If 
r'  contains  more  errors  than  the  code  can  correct,  then  the 
output  will  not  be  y  in  general.  Having  to  compute  a  higher- 
dimension  matrix-vector  product  is  traded  for  the  ability  to 
correct  errors. 

Our  work  to  date  has  involved  comparing  the  performance 
of  OMVMs  employing  error-correction  codes  to  OMVMs  not 
employing  the  codes  through  the  use  of  computer  simulations. 
The  results  of  this  research  have  been  published  in  [3]  and 
[4].  To  allow  for  a  mere  quantitative  assessment  of  the 
ability  of  the  codes  to  enhance  the  accuracy  of  an  OMVM,  a 
new  model  of  an  OMVM  that  incorporates  the  dynamic  ranges  of 
the  devices  has  been  developed.  Using  this  new  model, 
various  computer  simulations  have  been  performed  that  give  a 
quantitative  value  for  the  accuracy  increase.  Figure  2  shows 
one  simulation  result.  The  graph  was  obtained  by  using  the 
single-error-correcting  code  presented  in  [1].  For  a  signal- 
to-noise  ratio  (SNRy  )  of  approximately  8,  i.e. 
approximately  18  db.,  in  the  output  of  the  OMVM,  an  accuracy 
increase  of  approximately  .3  bits  can  be  obtained.  For  this 


ACCURACY  INCREASE  (BITS) 


STANDARD  DEVIATION  OF  NOISE  X  1 0^ 
SNRy  ->  19  9-5  6.4  4.8  3.8  3.2 


Figure  2.  Accuracy  increase  versus  signal-to-noise  ratio  in  output 


simulation,  a  strict  error  measure  was  ''mpioyed.  An  error  is 
said  to  have  occurred  if  any  component  of  the  product  vector 
is  incorrect.  By  using  a  less  stringent  error  measure,  such 
as  mean  square  error,  the  accuracy  increase  should  be 
greater.  Codes  capable  of  correcting  more  errors  should 
provide  an  even  greater  improvement.  Currently,  research 
continues  in  this  area. 

Various  techniques  have  been  introduced  over  the  past 
decade  to  achieve  high  accuracy  computation  with  analog 
optical  matrix- vector  multipliers  (OMVMs) .  The  two  most 
discussed  are  digital  multiplication  by  analog  convolution 
(DMAC)  and  the  residue  number  system.  Recently,  a  new 
technique,  digital  partitioning,  has  been  proposed  [5]  .  The 
basic  idea  behind  digital  partitioning  is  to  divide  a  large 
dynamic  range  matrix-vector  product  into  a  series  of  smaller 
dynamic  range  subproducts.  A  weighted  sum  of  these 
subproducts  will  give  the  large  dynamic  range  final  matrix- 
vector  product.  One  potential  drawback  of  digital 
partitioning  is  the  effect  of  small  errors  in  the  individual 
subproducts.  The  errors  can  occur  due  to  noise  in  the 
analog  system  and  also  due  to  component  failures.  A  small 
error  that  occurs  in  one  of  the  subproducts  that  is  then 
weighted  by  a  large  factor,  will  result  in  a  larger  error  in 
the  final  matrix-vector  product.  This  effect  can 
significantly  degrade  the  performance  of  the  digital 
partitioning  technique.  An  algorithmic  technique  that  can 


reduce  this  error  amplification  effect  is  the  application  of 
error-correction  codes. 

Error-correction  codes  reduce  the  error  ampl i f ication 
effect  of  the  digital  partitioning  technique  by  correcting 
the  small  errors  in  the  subproducts  and  thus  not  allowing 
them  to  be  weighted  by  a  large  weighting  factor  and 
resulting  in  a  larger  error  in  the  final  result.  The  codes 
accomplish  this  error-correction  by  the  introduction  of 
redundant  information  into  each  matrix-vector  subproduct  as 
discussed  above.  The  resulting  redundant  information  in 
each  vector  subproduct  allows  the  correction  of  some  errors. 
The  number  and  magnitude  of  errors  that  can  be  corrected  is 
dependent  upon  the  amount  of  redundant  information 
introduced.  We  are  currently  investigating  the  application 
of  error-correction  codes  to  combating  the  error 
amplification  effect  of  the  digital  partitioning  algorithm. 
Various  computer  simulations  are  being  perform.ed  that  show 
the  reduction  in  the  error  magnitude.  Figure  3  shows  the 
result  of  one  of  our  computer  simulations.  This  graph  is  a 
plot  of  the  root-mean-square  (rms)  error  divided  by 
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Signal-to-noise  ratio  in  detectors 

Figure  3.  Graph  of  rms  error  versus  signal  to  noise  ratio  in 
detector  for  the  cases  of  no  error-correcting  code 
and  a  single  error-correcting  multilevel  Hamming 
code . 


the  length  of  the  output  vector  versus  the  signal  to  noise 
ratio  in  the  detector.  One  curve  shows  the  performance  when 
no  error-correction  code  is  being  used  while  the  other  curve 
shows  the  performance  with  a  single-error-correcting 
multilevel  Hamming  code.  There  is  clearly  a  region  between 
very  low  and  very  high  noise  levels  in  which  the  error- 
correction  code  can  improve  performance. 


In  addition  to  the  research  mentioned  above,  we  are 
investigating  the  effects  of  scaling  the  codes  to  higher 
dimension  matr ix-vector  products.  Since  the  benefits  of  the 
codes  are  gained  at  the  expense  of  more  computations,  an 
analysis  of  the  additional  computations  required  is  being 
made.  Codes  more  suitable  to  OMVMs  are  being  developed  also. 
Some  codes  currently  used  have  the  drawbacks  of  requiring 
more  levels  to  represent  the  redundant  information  than  the 
actual  data  [1]  or  the  codes  require  that  a  residue  matrix- 
vector  multiplication  be  performed.  Optical  implementations 
of  the  encoders  and  decoders  are  also  being  investigated. 


2.  Statistical  Analv&is....and  Modeling  of  Optical  RrQceaiSQrs 


The  two  main  objectives  of  this  phase  of  our  research 
have  been  (1)  to  broaden  the  range  of  validity  of  our 
earlier  statistical  analysis  and  modeling,  and  (2)  to  apply 
the  general  analysis  to  specific  cases  of  interest  in  an 
effort  to  gain  a  more  pragmatic  view  of  the  system 
performance . 

Our  earlier  work  involved  a  general  statistical 
analysis  of  the  Stanford  optical  matrix-vector  multiplier 
[6]  [see  Figure  4]  under  certain  simplifying  assumptions, 
establishing  the  probability  density  and  mutual  coherence 
functions  of  the  detector  field  intensity  and  the  detector 
output  current  in  terms  of  the  statistical  characteristics 
of  the  system  components  [7]  .  Specifically,  we  considered 
the  effects  of  diffraction  leading  to  an  interchannel 
crosstalk,  carefully  analyzed  the  effects  of  the 
photodetection  process  [see  Figure  5]  on  the  signal 
statistics,  and  incorporated  various  postdetection  noise 
processes  [see  Figure  6]  into  our  final  results  [8,  9]  .  The 
expressions  produced  by  this  effort  proved  to  be  relatively 
complicated  and  generally  intractable  even  under  the  most 
idealized  conditions  [5,  61. 

One  of  our  goals  over  the  past  year  has  been  to  extend 
our  analysis  to  account  for  more  realistic  situations.  In 


Figure  4  Bl<x:k  diagram  of  the  Stanford  optical  matrix-vector  multiplier 


Figure  5 


Block  diagram  of  the  photodelcction  process 


F’igure  6  Block  diagram  of  the  post-detection  processing 


Figure  7 


Channel  geometry  for  crosstalk  analysis 


particular,  we  have  attempted  a  more  complete  optical 
analysis  of  the  system  by  bringing  in  the  crosstalk  and 

background  effects,  which  were  previously  mildly  studied  or 
altogether  neglected  (12,  13,  14].  We  used  the  more 
realistic  Rayleigh-Sommerfeld  diffraction  formula,  rather 
than  the  Fresnel  and  Fraunhofer  approximations,  to  be  able 
to  come  up  with  a  faithful  statistical  representation  for 
the  crosstalk  [see  Figure  7].  Background  light  fields,  on 
the  other  hand,  were  modeled  as  zero-mean  and  uncorrelated 
Gaussian  processes,  as  is  commonly  done  in  astronomical 

optics.  In  the  photodetection  and  postdetection  stages,  we 
extended  our  analysis  to  include  the  possibility  of 
(avalanche)  gain  noise,  and  also  brought  in  other  noise 
sources,  such  as  the  1/f  noise  and  the  generation- 

recombination  noise,  to  complete  our  work  [15]. 

Another  goal  in  our  investigation  was  to  put  our 
results  in  a  more  pragmatic  perspective  by  exploring  the 
practical  implications  of  these  theoretical  results. 

Specifically,  we  considered  classical  sources  such  as  lasers 
and  LEDs  for  which  the  device  statistics  have  been  well 
established.  Since  statistical  models  for  SLMs  are  currently 
unavailable,  we  considered  theoretical  models  such  as 
Gaussian  and  Bernoulli  for  the  intensity  transmittances  of 
these  devices.  Deterministic  modulation  effects  in  the 
sources  and  SLMs  can  easily  be  accounted  for  if  desired,  and 


were  neglected  here  for  brevity.  Among  the  photodetector 
filter  models  considered  were  ideal  ones  such  as  rectangular 
and  triangular  as  well  as  more  realistic  ones  such  as 
exponential  and  logarithmic  [16] . 

The  analytical  intractability  of  the  results  has 
prompted  us  to  resort  to  computer  simulations.  Currently, 
efforts  are  under  way  to  model  the  statistical 
characteristics  of  the  system  via  the  Monte-Carlo  technique 
to  produce  results  that  support,  and  supplant,  analytic 
solutions,  and  thus  offer  insights  into  the  practical 
significance  of  the  different  noise  sources  as  well  as  the 
practical  forms  of  the  output  signal  and  noise  statistics. 
We  are  also  hoping  to  gain  some  practical  insight  from  the 
results  of  the  experimental  phase  of  the  investigation,  (see 
Section  3)  which,  combined  with  the  simulation  results, 
should  enable  us  to  apply  promising  detection  and  estimation 
techniques  and  assess  the  practical  magnitude  of  the 
achievable  accuracy  enhancement  [17] . 

Ultimately,  we  aim  to  apply  this  general  strategy  of 
analysis  and  modeling  to  alternate  systems  and 


architectures . 


3 .  Device  Characterization  for  Analog  Optical  Computing 


This  phase  of  the  project  is  concerned  with 
characterizing  the  performance  of  coitjnon  devices, 
particularly  spatial  light  modulators  (SLMs),  which  are  used 
in  analog  optical  computing  systems.  The  performance  issues 
of  interest  are  those  which  impact  the  accuracy  of  the 
device  in  representing  quantized  numeric  states.  To  this 
end,  the  temporal  characteristics  have  been  of  primary 
interest.  Spatial  attributes  such  as  signal  uniformity  and 
channel  crosstalk  will  take  on  greater  significance  when  the 
temporal  characterization  is  complete. 

The  device  characterization  phase  has  been  subdivided 
into  an  experimental  and  an  analytical  effort.  The 
experimental  effort,  which  to  date  has  engaged  the  majority 
of  our  time,  is  concerned  with  developing  a  low  noise  test 
facility,  instrumentation,  and  data  acquisition  and  analysis 
tools.  The  analytical  effort,  which  is  currently  being 
given  more  attention,  is  concerned  with  developing  physical 
and  empirical  models  of  the  SLMs  so  that  device  behavior  can 
be  understood  and  meaningful  experiments  can  be  designed. 

The  overall  goal  of  the  device  characterization  is  to 
develop  a  complex  field  description,  through  physical  models 
and  experiment,  to  describe  the  accuracy  performance  of 
common  SLMs  used  for  the  representation  of  quantized  numeric 
information.  These  physical  models  and  experiments  can  then 
be  used  to  create  a  statistical  description  of  the  systems 


performances.  These  statistical  results  will  then  be 
readily  usable  in  Sections  1  and  2  of  this  project. 
Although  the  SLM  is  the  device  of  primary  interest  for 
modeling,  analysis  will  be  performed  on  sources  and 
detectors  as  well. 

We  shall  discuss  progress  on  each  of  the  experimental 
and  analytical  efforts  in  turn. 

3.1  Experimental  Aspect  of  Device  Characterization 

The  goal  of  the  experimental  effort  is  to  gather  high 
quality  data  (intensity  and  phase)  on  the  temporal  and 
spatial  device  characteristics  which  impact  the  accuracy  of 
represented  quantized  numeric  data  for  devices  commonly  used 
in  optical  computing.  Temporal  stability  of  a  given  signal 
state  and  the  response  times  between  signal  levels  are  of 
current  experimental  interest.  Although  the  basic 
performance  attributes  of  a  device  or  family  of  devices  are 
of  most  interest,  signal  artifacts  resulting  from  the  drive 
electronics  packaged  with  the  devices  must  also  be 
characterized.  In  general  the  devices  of  experimental 


interest  include: 


_ Modulators _ 

Liquid  crystal  (LC)  (both  nematic 

and  ferroelectric) _ 

Acousto  -optic  (AO) _ 

Magneto-optic  (MO) _ 

_ Sources _ 

Helium-Neon  (HeNe) _ 

Argon  Ion  (Ar+) _ 

Light  emitting  diodes  (LEDs) _ 

Laser  diodes _ 

_ Detectors _ 

PIN  photo  diodes _ 

Avalanche  photo  diodes _ 


Table  I.  General  device  classes  of  experimental  interest. 

3.1,1  Facilities 

To  realize  the  experimental  goals,  a  low  noise  test 
facility  along  with  the  associated  instrumentation  and  data 
acquisition/analysis  tools  has  been  built.  The  facility 
consists  of  a  20'xl0'  Faraday  screen  room  (Figure  8)  which 
provides  a  low  radio  frequency  interference  (RFI) 
environment  within  which  sensitive  measurements  can  be  made. 
Over  the  past  year  shield  integrity  has  been  improved. 


Additionally,  component  shielding  of  power  supplies  and 
detector/amplifier  circuits  within  the  screen  room  has  been 
undertaken  to  minimize  the  noise  resulting  from  such 
internal  RFI  noise  sources  as  laser  power  supplies. 

Ai*  Beam  Shielded 


Figure  8.  The  RFI  screen  room  facility. 

Along  with  quieting  the  RFI  environment,  we  have  also 
fabricated  a  second  generation  of  low  noise  photo-detector 
circuits.  These  detector  circuits  include  a  prototype  of  an 
electronic  noise  canceller  that  has  the  potential  to  reach 
the  shot  noise  limit[18].  The  design  and  fabrication  of 
these  detector  circuits  has  been  motivated  by  the  severity 
of  laser  noise  revealed  by  our  early  experiments.  We  will 
discuss  some  of  the  measurement  techniques  and  rationale  for 
these  detector  circuits  in  the  section  on  experimental 
methods  and  initial  results. 


3.1.2  Data  Acquisition  and  Analysis  Tools 

Within  the  screen  room  facility  we  have  a  33Mhz,  386- 
based  PC  with  8Mbs  of  RAM.  For  data  acquisition  this 
machine  contains  an  Analogic  Fast  14-4,  IMhz  ADC  with  4Mb  of 
capture  memory.  Over  the  past  year  we  have  developed  the 
Optical  Analysis  Simulation  Interactive  System  (OASIS), 
which  is  a  Microsoft  (MS)  Windows  3.1  application.  OASIS 
has  been  built  to  be  the  main  data  acquisition,  analysis, 


General  Features  of  OASIS 

Support  for  20  channels  of  10,240  samples  per 
channel. (the  samples  per  channel  are  limited  by 
the  number  of  active  channels  and  the  amount  of 
available  RAM  on  the  host  system) 

File  support  for  reading  data  files . (support  for 
saving  files  will  be  added  shortly) 

Statistical  Analysis 

Compilation  of  histograms  with  bin  resolutions  from 
512  to  10240  over  dynamic  range  of  data  for  all 
channels.  Compilation  of  several  moments. 

Hypothesis  testing  between  channels  and  with 
synthetic  distributions  based  on  the  Chi  Squares 
and  Kolmogorov-Smirnov  tests. 

Auto  and  Cross  Correlation  and  Power  Spectral 

Density  analysis  between  channels. 

Frequency  Domain  Analysis 

Forward  transforms  with  several  time  domain 
windowing  functions. 

A  wide  range  of  frequency  domain  filters  (magnitude 
and  phase)  which  can  be  constructed  from  primitive 
windowing  functions.  A  few  simple  deconvolution 
techniques 

Time  Domain  Operations 

A  Synthetic  Data  Generator  (SDG)  from  which 
complicated  signals  can  be  created  from  primitive 
functions,  including  several  noise  models. 

Time  domain  modeling  via  least  squares  and  cubic 
spline  fits  to  the  channel  data. 

Table  II.  Functionality  of  the  OASIS  software  system. 


filtering  and  modeling  tool  for  our  experiments.  The 
current  functionalities  of  OASIS  are  listed  in  Table  II. 

The  OASIS  system  is  graphically  interactive  and  has 
been  built  completely  in-house  using  the  MS  C7 . 0  compiler. 
This  in-house  development  al’ows  OASIS  to  be  augmented  in 
any  manner  as  the  experiments  and  analysis  progress . 
Examples  of  the  time  domain  data  window,  histogram  window, 
correlation  -  power  spectral  density  window,  and  frequency 
domain  window  are  provided  in  Figures  9  through  12. 

At  this  time  the  data  acquisition  is  performed  by  a  DOS 
application  {also  built  in-house)  and  saved  to  a  data  file. 
The  short  term  necessity  for  this  external  program  has  been 
the  difficulty  in  writing  an  MS  Windows  device  driver  for 
the  Analogic  hardware.  With  the  introduction  of  the  MS 
Windows  NT  operating  system  in  early  1993  this  task  will 
become  tractable  and  the  acquisition  function  will  be 
written  into  OASIS. 
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Figure  9. 

The  OASIS  main  data  window  showing  1000  points  of 
the  8192  points  for  two  active  data  channels.  This 
data  is  of  a  noisy  DC  power  supply (the  smoother  trace) 
driving  an  LED  and  the  detected  output  of  the  LED (the 
noisier  trace) .  For  this  case  the  power  supply  ripple 


is  the  dominant  signal.  Notice:  the  phase  difference 
between  the  DC  signal  and  the  LED  output  is  due  to 
unsynchronized  data  captures  and  is  not  a  real 
attribute  of  the  LED. 


Figure  30.  The  OASIS  statistics  window  showing  a  1024 
bin  histogram  of  the  two  data  channels. 


Figure  11.  The  OASIS  correlation  and  PSD  window 

showing  the  cross  correlation  (CCF)  of 
the  DC  power  supply  and  the  LED  output 
(top)  and  the  cross  spectral  density 
(bottom) .  Since  the  data  consists  of 
8192  points  the  CCF  and  PSD  are  also 
8192  points.  The  displayed  plot  windows 
have  been  centered  with  the  zero  of  the 
CCF  and  DC  component  of  the  PSD  in  the 
center.  The  fundamental  frequency  is  at 
120Hz. 


Figure  12.  The  OASIS  frequency  domain  window  showing 
magnitude  of  the  transform  for  the  DC  power 
supply  and  the  LED  output  (top)  and  the  phase 
(bottom) .  Since  the  data  consists  of  8192 
points  the  transforms  are  also  8192  points. 
The  displayed  plot  windows  have  been  centered 
with  the  DC  component  in  the  center.  The 
fundamental  frequency  is  at  120Hz 


3.1.3  Experimental  Methods  and  Current  Results 


To  date  the  majority  of  experiments  have  been  conducted  to 
isolate  problems  in  the  facilities,  tools,  and  equipment.  A 
primary  concern  borne  out  by  these  experiments  has  been  the 
significant  level  of  noise  identified  on  both  the  HeNe  and 
Ar+  lasers.  This  noise  results  from  many  deterministic  and 
random  factors  ranging  from  power  supply  ripple  and  acoustic 
vibrations  to  longitudinal  mode  and  transverse  mode  lopping. 
We  have  systematically  proceeded  to  isolate  and  reduce  these 
effects  by  several  methods,  but  have  not  as  yet  been 


completely  successful  lu  the  exclusion  of  this  noise. 
Although  these  attempts  to  quiet  the  sources  is  continuing 
it  is  unlikely  that  this  noise  will  be  elim.inated  as  an 
experimental  concern.  In  conjunction  with  these  efforts  we 
have  also  pursued  methods  of  post-process  removal  of  this 
source  noise  from  the  experimental  results.  One  method 
currently  used  relies  on  concurrent  sampling  of  two 
detectors,  one  detecting  the  unperturbed  laser  beam  and  the 
other  the  output  of  the  SLM  and,  using  the  correlation  and 
frequency  domain  techniques  of  the  OASIS  system,  post¬ 
process  to  remove  the  laser  effects  [19].  To  date  this  has 
been  qualitatively  sufficient.  New  post-processing 
techniques  are  currently  under  development.  Another 
approach  to  removing  the  laser  noise  has  been  to  use  an 
electronic  noise  canceller  which  simulta.neously  detects  the 
source  and  SLM  beams  and  performs  an  electronic  analog 
subtraction [ 1 8 ] .  This  module  has  been  fabricated  and  is 
currently  under  test  and  refinement. 

In  addition  to  noise  on  the  laser  sources,  significant 
ripple  has  been  identified  on  the  power  supplies  driving 
the  detector  and  LED  circuits  (as  seen  in  the  Figures  9-12). 
Several  methods  have  been  satisfactorily  employed  to  reduce 
this  noise  to  acceptable  levels.  The  data  depicted  in  the 
figures  represents  the  state  of  the  system  before  these 
techniques  were  implemented. 


As  the  experimental  facility  has  evolved  we  have  been 
actively  testing  two  commercially  available  SLMs:  the  Hughes 
4050  Liquid  Crystal  Light  Valve  (LCLV) ,  and  the  Casio  TV- 
1400  Liquid  Crystal  Television  (LCTV) .  Early  results  show  a 
significant  level  of  temporal  fluctuation  on  both 
devices[19].  These  fluctuations  are  predominantly  due  to 
their  drive  electronics. 

To  fully  realize  the  limitations  of  an  LC-based  device 
for  the  application  of  representing  quantized  numeric 
information  we  must  have  access  to  the  LC  cell  without  the 
limitations  of  the  drive  electronics.  To  realize  this  we 
are  currently  planning  to  acquire  test  LC  cells  which  have 
no  prepackaged  drive  electronics.  These  cells  will  be 
fabricated  using  standard  orientations  of  the  LC  material 
found  in  commercial  devices  (i.e.  45,  90,  and  270  degree 
twists) .  With  these  cells  we  will  then  be  in  a  position  to 
supply  a  drive  signal  of  our  choosing  and  to  control  and 
meter  all  aspects  of  the  device.  In  this  way  we  can  develop 
quantifiable  bounds  on  the  performance  of  the  LC  cell 
without  being  limited  by  non-optimal  drive  electronics. 
These  closer-to-ideal  devices  will  play  prominent  roles  in 
the  analytical  phase  of  this  project  since  all  device  and 
material  characteristics  will  be  known.  With  appropriate 
ittodels  and  performance  bounds  established  we  can  then  go  to 
the  commercial  devices,  with  their  unknown  parameters  and 


non-optimal  drive  electronics,  and  perform  meaningful 
experiments  and  analysis. 

3.2  Analytical  Aspects  of  Device  Characterization 

Although  the  bulk  of  our  effort  has  been  focused  on 
getting  the  experimental  tools  in  place  we  have  recently 
turned  more  effort  toward  the  analytical  aspects.  Here,  the 
basic  understanding  of  the  physics  of  the  device  under  test 
is  fundamental  to  the  design  of  the  experiments. 

Our  first  step  has  been  to  define  "noise"  and 
"accuracy"  from  the  device  perspective.  The  crux  of  the 
problem  lies  in  the  fact  that  we  want  to  represent  quantized 
numeric  information  on  basically  linear  devices. 
Understanding  the  behavior  of  the  device  at  specific  signal 
levels  now  becomes  the  mechanism  for  determining  the 
achievable  dynamic  range  and  bit-resolution.  Knowing  the 
physics  of  the  device  not  only  allows  one  to  understand  it's 
behavior  but  it  also  provides  insight  into  what  can  be 
changed  to  increase  it's  performance.  One  of  the  analytical 
goals  of  this  project  is  to  provide  this  analysis  and  to 
produce  means  by  which  different  device  construction  and/or 
operation  can  achieve  better  accuracy  performance. 

The  analytical  phase  will  develop  along  two 
complimentary  paths.  The  first  is  an  investigation  into  the 
device  physics.  To  this  date  the  focus  has  been  on  nematic 


liquid  crystals  (NLCs)  (both  twisted  and  non- twisted) ,  and 
ferroelectric  liquid  crystals  (FLCs) .  It  is  envisioned  that 
a  dynamical  model  for  the  LC  materials  related  to  accuracy 
will  be  developed.  This  model  will  be  used  to  guide 
experiments  and  to  predict  conditions  in  which  the  LC-based 
SLM  will  perform  better.  The  second  avenue  is  along  a  purely 
empirical  investigation.  The  goals  of  this  will  be  to  1) 
establish  the  predominant  dynamical  features  of  the  SLM 
without  being  rigorously  exact,  and  2)  utilizing  system 
identification  techniques,  develop  the  best  communications 
system  model  for  the  "black  box"  SLM.  Both  of  these  results 
will  directly  augment  the  device  physics  models  and  drive 
the  experimental  investigation. 

Other  classes  of  devices  such  as  AO  and  MO  modulators 
will  be  investigated  in  a  similar  manner  as  the  project 


progresses . 
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1.  Development  and  simulation  of  a  system  model 
incorporating  device  dynamic  range  for  better  quantitative 
assessment  of  error-correction  code  performance. 

2.  Extension  of  earlier  statistical  models  to  include 
crosstalk,  background,  avalanche  gain,  flicker  and 
generation-recombination  noise  effects. 

3.  Analytical  and/or  numerical  evaluations  of  the  output 
signal  statistics  for  a  number  of  practically  interesting 
device  combinations. 

4.  Construction  of  the  Optical  Analysis  Simulation 
Interactive  System^  (OASIS)  software  for  the  acquisition, 
analysis  and  manipulation  of  experimental  data. 

5.  Identification  of  major  noise  sources  of  experimental 


concern  using  OASIS. 


